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ABSTRACT 



We test a nonlinear force-free field (NLFFF) optimization code in spheri- 
cal geometry using an analytical solution from Low and Lou. Several tests are 
run, ranging from idealized cases where exact vector field data are provided on 
all boundaries, to cases where noisy vector data are provided on only the lower 
boundary (approximating the solar problem). Analytical tests also show that the 
NLFFF code in the spherical geometry performs better than that in the Carte- 
sian one when the field of view of the bottom boundary is large, say, 20° x 20°. 
Additionally, We apply the NLFFF model to an active region observed by the 
Helioseismic and Magnetic Imager (HMl) on board the Solar Dynamics Obser- 
vatory (SDO) both before and after an M8.7 flare. For each observation time, 
we initialize the models using potential field source surface (PFSS) extrapola- 
tions based on either a synoptic chart or a fiux-dispersal model, and compare the 
resulting NLFFF models. The results show that NLFFF extrapolations using 
the flux- dispersal model as the boundary condition have slightly lower, therefore 
better, force-free and divergence-free metrics, and contain larger free magnetic 
energy. By comparing the extrapolated magnetic field lines with the extreme 
ultraviolet (EUV) observations by the Atmospheric Imaging Assembly (AIA) on 
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board SDO, we find that the NLFFF performs better than the PFSS not only for 
the core field of the flare productive region, but also for large EUV loops higher 
than 50 Mm. 

Subject headings: Sun: activity — Sun: corona — Sun: magnetic topology — 
Sun: surface magnetism 



1. Introduction 



Magnetic field plays an important role in structuring the thermodynamic and hydrody- 
namic behaviors of the plasma in the solar atmosphere. Three dimensional magnetic fields, 
especially in the higher solar atmosphere, provide crucial information toward understand- 
ing various solar activities, such as filament eruptions, flares, and coronal mass ejections 
(CMEs). However, there are some difficulties in measuring the coronal magnetic field. The 
emissions and field strength are much weaker in the corona than that in the photosphere 
and chromosphere, while the temperature is much higher producing a wide line width that 
covers the splitting by the Zeeman effect. Additionally, the observed intensities are in- 
tegrations along the line-of-sight direction due to the optically thin nature of the coronal 
plasma, and thus bright photospheric and chromospheric features obscure the emission origi- 
nating in the corona. Although measurement of magr ietic fields in the chrornosphere and the 
coron a has considerably improved in recent decades ( 1 Judge! Il998l : iLin et al.ll2000l : iLiu fc Lin 
20081 ). further developments are needed before accurate data are routinely available. Since 
photospheric magnetic fields are widely observed by various ground-based and space-borne 
instruments, and since the magnetic fields are in a low (5 (the ratio between the gas pressure 
and the magnetic press ure) environ ment in the upper chromosphere, the transition region, 
and the lower corona f lGaryl l200ll ). a force-free field extrapolation technique provides an 
alternative way to infer the magnetic fields in the upper solar atmosphere. 



The force-free magnetic field is described by the equations as follows: 



V X B 
VB 



aB 

, 



(1) 
(2) 



where a is the torsion function that represents the proportionality between the electric cur- 
rent density and the magnetic field. Combining equations ([T]) and ([2]), one obtains B-V« = 0, 
which indicates that a is always constant along a given field line. Setting a = yields the 
current-free potential field. For constant a over space, equations ([1]) and ([2]) describe the 
linear force-free field. In the general nonlinear force-free field (NLFFF) that is closer to re- 
alistic coronal magnetic fields, a varies from one field line to another. Solutions to the fully 
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nonlinear problem can be calculated using numerical algorithms such as the Grad-Rubin, up- 
ward integration, MHD relaxatio n (magnetofriction al) , optimization, and boundary element 
(Green's function like) methods (jWiegelmannll2008l and references th erein). A limited class 
of axisymmetric semi-analytical NLFFFs also exist (ILow fc Loulll990[ ). and these are usually 
adopted as t he standard test mod el of the above mentioned numerical algorithms, as has 
been done in lSchrijver et al.l (120061 ) . Many authors have shown that the force- free function a 
appe aring in Equation ([1]) changes in active regions (e.g.. lRegnier et al.ll2002l : ISchrijver et al. 



20051 ). thus the NLFFF model is necessary to construct the magnetic configuration for active 
regions, especially those with eruptive activities. 

NLFFF models have been adopted to st udy various rnagneti c field structures and prop- 
erties in the solar atmosphere . For instance. ICanou et al.l (|2009[ ) found a twisted flux rope 
before eruptive activities and ICanou fc Amaril (120101 ) obtained a flux rope as the magnetic 
structure of a filam ent using the Grad-Rubin method. With the MHD relaxation method, 
Bobra et al.l (120081 ) and ISu et al.l ( 120 111 ) constructed non-potential magnetic field models for 
solar active regions and a fiare/CME event, respectively. The MHD relaxation method is 
implemented in spherical geometry. A flux rope is usually present in this model si nce the flux 



rope is inserted into a potential field beforehand. Using the optimization method, iGuo et al. 



( l2010al Jbl) found that a flux rope and magnetic arcades coexisted in a solar filament and 
they fu rther computed the twi st number of the flux rope and studied the eruption mech- 



anism. 



Jing et al.l (120091, l2012l ) studied the evolution of free magnetic energy and relative 



magnetic helicity in activity productive active regions with the same optimization method. 
It is worth mentioning some extrapolation methods without the force-free field assumption, 
for example, the non-force-free extr apolation of co r onal m agnetic field based on the principle 
of minimum dissipation rate (e.g.. iHu et al.ll2008l. 120101) and the magnetohydrostatic mod- 



eling of the corona (e.g., IWiegelmann fc Inhesterll2003t IWiegelmann et al.l 120071 : iRuan et al. 



20081 ). Those methods should be more applicable when the force-free state is no longer valid. 



Most of the NLFFF procedures are implemented in the Cartesian coordinates. These 
procedures are not well suited for larger domains, since the spherical nature of the solar 
surface cannot be neglected when the field of view is too large. However, in a critical test 
of NLFFF codes, iDeRosa et al.l (120091 ) concluded that the field of view should be as large 
as possible to include more information of magnetic field line connections. In addition to 
the above reason, when we apply the optimization method for the NLFFF extrapolation to 
any realistically observed active region, an apriori requirement is that the magnetic field on 
the boundary is isolated ( Wiegelmann et al.ll2006l ). However, this requirement may not be 
satisfied in a realistic observation due to the limitation of the field of view. Therefore, it 
appears prudent to implement a NLFFF procedure in spherical geometry for use when large- 
scale boundary data are available, such as from the Helioseismic and Magnetic Imager (HMI; 
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Scherrer et al.l l2012l : ISchou et al.l |2012| ) on board the Solar Dynamic Obseratory (SDO). 
Furthermore, studies of magnetic field structures of large scale solar phenomena, such as 
trans-equatorial magnetic loops, large filaments, and CMEs, are also likely to benefit from 
NLFFF procedures in the spherical geometry. 

A NLFFF proce dure in the spheric al geometry based on the optimizatio n method has 



been i mple mented bvlWiegelmannI (120071 ). which has been further developed by lTadesse et al. 



(120091 ) and lTadesse et al.l (120111 ). Another NLFFF procedure in the spherical geometry based 
on the same method has been implemented by J. McTiernan in the FORTRAN 90 lan- 
guage and released in the "nlfff" package available through Solar Software (SSW^^^Afthough 
adesse et al.l (120091 ) have already tested the NLFFF code implemented by IWiegelmann 



(120071 ) with the Low and Lou analytical solution, we test again J. McTiernan's version of 



the NLFFF code in this paper carefully. It is not only meaningful for the code itself before 
we apply it to real observations, but also an independent validation for the NLFFF method 
in spherical geometry. 

We then apply the NLFFF procedure to a group of active regions observed on 2012 
January 23 by SDO /HMl. An M8.7 class fiare occurred in NOAA 11402, one among the 
active region group. With this study, we first compare the extrapolated (both potential 
and NLFFF) magnetic loops wi th extreme ultravio let (EUV) observations by the Atmo- 
spheric Imaging Assembly (AIA; iLemen et al.ll2012l ) on board SDO. This comparison indi- 
cates whether the NLFFF model reconstructs the magnetic configuration better than the 
potential field model. Additionally, this comparison can be used to evaluate how well the 
model field lines approximate the observed coronal loops, thereby enabling a determination 
of whether the model is consistent with observations. Furthermore, the spherical modeling 
technique enables the study of the non-potentiality of large magnetic loops and the magnetic 
surroundings of an active region, especially in a place where many active regions interact 
with each other. Finally, such a study can show how a fiare, together with the associated 
filament eruption and CME, changes the magnetic configuration permanently both on the 
photosphere and in the corona. 

This paper is organized as follows. We first describe the procedure implementation in 
Section [21 The testing of the procedure using the Low and Lou analytical solution, and 
comparisons between the NLFFF extrapolations in the Cartesian and spherical coordinates 
are presented in Section [31 Then, we apply the NLFFF extrapolation in spherical geometry 
to SDO /B.MI observations and compare the NLFFF models with EUV loops observed by 
SDO /AIA in Section [H A summary and discussions are finally presented in Section [SI 
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2. Optimization Procedure in Spherical Geometry 



Wheatland et al.l f l2000l ) proposed an optimization method to reconstruct the NLFFF 
by minimizing an objective functional that combines Lorentz forces and the divergence of 
the magnetic field. If the functional is minimized to zero, equations ([1]) and ([2]) are satisfied 
si multaneously The op timi zation procedure in th e spherical geometry has been implemented 



by IWiegelmannI (120071 ) and iTadesse et al.l (120091 ) . where the objective functional is defined 



as 



L^= [ uj{r,e,^) \b-^ I (V X B) X B |2 + I V ■ B l^lr^sinMrdM^ , (3) 
Jv ^ 

where w{r, 6, </>) is a weighting function and V is a. wedge-shaped volume for numerical 

computation. The weighting function is adopted to minimize the effects of the lateral and 

top boundaries, where the data are not available in practical observations. The wedge-shaped 

volume is defined by six boundaries with r e [r^in, rmax], 6 e [6'min, ^'max], and € [0min, 0max]- 

The whole computational box is resolved by the grids of Ur x ng x n^. 

For test cases with ideal boundary conditions, the vector field on all six boundaries can 
be specified by analytical solutions, which requires equal weights for all grid points in V with 
u{r, 6,(j)) = 1 everywhere. However, in more realistic cases, the vector field is known only on 
the lower boundary. In these instances the boundary conditions on the remaining four sides 
and top boundaries are typically specified by potential fields based on extrapolations of the 
photospheric line-of-sight field. Because the potential field on the five boundaries does not 
match the NLFFF assumption exactly, the objective functional L^^ cannot decrease to zero. 
For these non-ideal cases we introduce a buffer region for the lateral and top boundaries, 
where the weighting function decreases from 1 to with a cosine profile, from the inner 
volume V to the edges, in the buffer region with ribuf grid points. Thus, the inner physical 
domain V with weights of 1 spans (n^ — nbuf) x {ng — 2nbuf) x {nif, — 2nbuf) grid cells. The 
weighting function above enables departures of the field from the force-free and solenoidal 
state in the buffer region, which allows the magnetic field in the inner physical domain V 
to evolve to a more force-free state than the code without a weighting function does. This is 
because that those field lines with one of their footpoints on the lateral and top boundaries 
are half free, i.e., being allowed to change during the optimization procedure. The final 
state of the force-free field depends on the initial condition and the iteration history of the 
optimization procedure. 

In this paper, we experiment with the optimization procedure in the Solar Software 
(SSW), which is implemented both in the Cartesian and spherical geometries and includes 
the weighting function, to reconstruct large scale coronal magnetic fields. This version of the 
optimization procedure that is fulfilled in the Cartesian coordinate system without weighting 
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function has been tested in ISchrijver et al.l (120061 ) and iMetcalf et al.l ( 120081 ) . The author of 
the code is J. McTiernan. In the version that we adopt here, the procedure includes the 
weighting function as an input file. 



3. Testing the Optimization Procedure with the Low and Lou Solution 

In this section, we test the optimization procedure as described in Section [2] using the 



Low fc Loul (Il990l ) solution. The subsections are arranged as follows. We first describe the 
settings of the Low and Lou solution that is used as the standard test model in Section 13.11 
Secondly, the potential field source surface (PFSS) models that are used as initial conditions 
for the optimization method are discussed in Section [3^ Thirdly, our testing results both for 
the ideal and noisy boundary conditions are presented in Section [3l3l Finally, we compare the 
extrapolation results derived by the NLFFF code in the Cartesian coordinates and spherical 
coordinates, respectively, in Section 13.41 



3.1. Low and Lou Solution 



The Low and Lou solutions describe a set of NLFFFs in axially symmetric geometries 
associated with point sources placed at the origin. Let us define the heliocentric coordinate 
system that is denoted by 0-xyz (also referred as the physical coordinate system), where O 
is the center of the Sun and xyz are the Cartesian components. If the origin of the local 
coordinate system, O'-XY for the Low & Lou solution is placed at (— /x;0,— Z^) and the 
local F-axis is parallel to the physical y-axis, the relation between the local coordinates and 
the physical ones is given by 



X = {x + Ix) cos ^ - {z + Iz) sin $ 
Y = y , 

Z = (x + Ix) sin $ + (z + cos $ 



(4) 



where $ is the angle counterclockwise measured from z-axis to Z-axis if viewed against y-axis. 
Al though the t r ansfo rmation is modified a little as shown in Equation |4] from Equation (13) 
of ILow fc Loul (119901) ■ the coordinate system is identical to that depicted in Figure 2 of 



Low k Loul (Il990h . 



We adopt the NLFFF solution with the eigenfunction of Pi i in the local coordinate 
system. Pi.i is a solution of a nonlinear ordinary differential e quatio n, corresponding to 
parameters n = 1 and m = 1 of Equation (5) of Low fc Lou ( 199ol ). The final Low & 



Lou solution in the physical volume with r G [l.Oi?©, 2.5i?0], 9 G [0.0°, 180.0°], and 
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[0.0°, 360.0°] is calculated with the parameters Ix = O.25-R0, Iz = 0, and $ = — tt/IO, and 
resolved by 30 x 60 x 120 grid points. All the three components Bj., Bq, and B^ have been 
multiplied with a normalization constant to ensure that the maximum of \Br\ equals 500 
G. Some sample field lines of the Low & Lou solution are shown in Figure [DJ^a), where the 
central view point is located at [6', 0] = [90°, 180°], i.e., an observer is observing along the 
X-axis. 



3.2. Potential Field Source Surface Model 



A potential (current-free) field obeys the equation V x B = 0, so that the magnetic 
field can be expressed as the gradient of a scalar potential , i.e. B = — V^l '- Sinc e V ■ B = 0, 
the scalar potential obeys the Laplace equation V^\E' = 0. ISchatten et al.l (119691 ) introduced 



a spherical source surface at a distance Rs from the Sun's center, where the field lines are 
radial, to simulate the effect of the solar wind on the magnetic field. Therefore, the field 
vectors are purely radial at Rg, which indicates that \l/ is a constant on the source surface. 
This constant can be selected as 0. Together with the photospheric boundary condition, 
which is traditionally provided by a map of the radial magnetic field, the Laplace equation 
can be solved in the spherical coordinate system (r, 6, (p), where 6 stands for colatitude. The 
solution in the domain Rq < r < R^ is 



^(r, ^, 0) = £ Yl i^Tr^ + 5™r-('+^)]Y,'^ 

l=Om=-l 



(5) 



where the coefficients A]^ and -B™ are determined by the boundary conditions at Rq and Rs. 
The spherical harmonic functions Yi"^{6, 0) = C/"P™ (cos ^)e*™''^, where are normalization 
constants based on / and m, and P™(cos 9) are associated Legendre functions. Please refer 
to the appendix in ISchrijver fc De Rosal ( 120031 ) for a detailed deduction of the coefficients 
and the final magnetic field. Further in formation about spher i cal h armo nic expansion of 



the La place equation can also be found in lAltschuler fc NewkirkI ( 119691 ) and lAltschuler et al. 



( 119771 ). Using the radial magnetic field observed at Rq and the source surface assumption 
at -Rs as the boundary condition, together with the spherical harmonic expansion in the 
domain Rq < r < Rs, we finally obtain the PFSS model. Throughout this paper, the 
harmonic coefficients of the PFSS model are computed by the "pfss" package available in 

ssw. 

For the purpose of initializing the test runs based on the Low and Lou solution, we 
construct PFSS fields using the radial component of the magnetic field from the Low and 
Lou solution at the lower boundary. These fields are shown in Figure [U^b). The principal 
order of the spherical harmonic series, /max? is set to 7. The solutions are smooth enough that 



- 8 - 



higher / values are not needed to resolve the functional form of the Low and Lou solutions 
at these radii. We compute the correlation coefficient between the computed radial field (6^) 
and the analytical radial field {Br) and the normalized error (defined as ^ \br — Br\/ ^ \ Br\) 
on the bottom boundary as a measurement of the agreement between them. The correlation 
coefficient and normalized error of br and Br are 1.00 and 0.05, respectively, which indicates 
that the radial components of the analytical magnetic field are well reconstructed by the 
PFSS model. 



3.3. Test Results with Ideal and Noisy Boundary Conditions 

In this section, we will perform numerical experiments with the optimization procedure 
written by J. McTiernan in spherical geometry, with the aim of recovering the Low and Lou 
analytical solution described in Section [3?T] using different kinds of boundary conditions. Each 
experiment is initialized with potential fields calculated using the PFSS model described in 
Section 13.21 The NLFFF extrapolations are performed in a wedged-shaped domain bounded 
by r G [I.ORq, 2.5Rq], 9 e [9.0°, 171.0°], and G [90.0°, 270.0°]. This domain is resolved by 
30 X 54 X 61 grid points. For the ideal boundary test cases, we arrange the following three 
different boundary conditions: 



Case 1: The analytical vector magnetic fields from the Pi i solution of iLow fc Loul (jl990[ ) 
are specified on all the six surfaces for the boundary condition. 

Case 2: The Low and Lou fields are specified only on the bottom layer, and the others are 
specified by the PFSS. No weighting function is used for lowering the boundary effect. 

Case 3: With similar boundary conditions used in Case 2, but a buffer zone with 6 grid 
points toward the top and lateral boundaries is adopted, where the weighting function 
decreases from 1 to with a cosine profile. The computation domain is the same size 
as Case 2, but with the buffer inside the Case 2 boundary. 

In order to measure the convergence degree of the computed magnetic field, we define 
the integral of the Lorentz force and the magnetic field divergence divided by the volume as 
follows: 

Lf = 1 / w(r,^,0) I (V X B) X B |2 r^sinMrdM^ , (6) 
^ Jv 

Ld = ^ I w(r, ^, 0) I V ■ B |2 sin MrdM0 . (7) 



V 



V 
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The magnetic field and the length units are Gauss and Mm, respectively. The two integrals 
should be for perfect force-free and divergence-free field. The objective functional defined 
in Equation divided by volume V is the summation of Lf and L^, i.e., = V{Lf + Ld). 



We use the figures of merit defined in ISchrijver et al.l ( l2006l ) and iMetcalf et al.l (120081 ) 



to quantify the degree of agreement between the analytical magnetic field B and the extrap- 
olated b. The vector correlation metric is defined as 



- 

the Cauchy-Schwartz metric: 

Ccs 

the normalized vector error metric: 



1 \ ^ Bj ■ bj 



En = EilJ^» ^^1 ^ (10) 



and the mean vector error metric: 

Em 



1 b,' — B,- ^ , 

_\ li_ (11) 



Here, Bi = |Bj| and bi = |bj|. The summations run over all the points i in the volume of 
interest, and M is the total number of points. If the two vector fields are identical, then 
Cvec = 1 and Cqs = 1; C'vcc = and Ccs = if Bj ± bj at each point. Unlike the vector 
correlation and Cauchy-Schwarz metrics. En = and E^ = if the two vector fields are 
perfectly matched with each other. For a consistent comparison with the other metrics, we 
list 1 — En and 1 — Ej^ in the following tables, so that all the metrics reach unity if the two 
vector fields are identical. Additionally, the magnetic energy in magnetic field b normalized 
to that in the analytical magnetic field B is defined as 

where e = 1 when the two fields are identical. 

All the aforementioned metrics for the Low and Lou solution, the PFSS model with 
^max = 7, and the three test cases with ideal boundary conditions are listed in Table [TJ The 
total energy contained in the corresponding magnetic field and that contained in the Low 
and Lou solution is also listed in the last column of Tabled! The PFSS model has lower Cvec 
and Ccs nietrics, but larger E'n and Em metrics compared with the Low and Lou solution. 
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The PFSS magnetic energy is only 82% of the reference one. However, using the potential 
field as the initial condition, the NLFFF with all the six boundary vector magnetic fields are 
provided (Case 1) is recovered to the Low and Lou solution with very high accuracy as shown 
in Table [U For a more practical case, in which only the bottom vector magnetic fields are 
known, while the others are specified by the PFSS model, and we do not use the weighting 
function (Case 2), the results show that Lf and are one order of magnitude larger than 
that of the Low and Lou solution. However, the C^ec, Cqs, -E'n, and Em metrics are better 
than that of the PFSS model. The magnetic energy in Case 2 is recovered to 92% of the 
energy in the Low and Lou solution. If a weighting function decreasing from 1 to with a 
cosine profile is adopted in the buffer zone (Case 3), Lf, L^, and the magnetic energy have 
similar values as that in Case 1. The C^ec, Cqs, -^'N, and Eu metrics are worse than that in 
Case 1 but better than that in Case 2. Therefore, a weighting function for those boundaries 
where data are missing improves the NLFFF extrapolation. 



Molodenskvl fll969l . Il974h . lAlvl (119891 ) . and lSakurail fll989f l pointed out that vector mag- 



netic fields on a closed surface that fully encloses any force-free domain have to satisfy the 
force-free and torque-free conditions. In practical observations, only data on the bottom 
boundary are available. Therefore, the force-free and torque-free conditions are required 
in a well isolated region in a force-free magnetic field. In an isolated region, all field lines 
originating from the bottom boundary fall back on it again. The isolation condition requires 
that the flux is in balance, an area with strong magnetic fields (active region) is surrounded 
by weak fields (quiet Sun), and no flux interconnects to other active regions. However, the 
flux balance is a necessary condition for the isolation of the magnetic field, but it is not a 
sufficient one, which implies that it is better to enlarge the field of view than to cut a smaller 
region in order to fulfill the isolation condition. On the photosphere the force-free and 
torque-fr ee condition s are not well fulfilled in an isolated region since the plasma /3 is close 
to unity (IGarvl 120011). T hus, the forced photospheric boundaries need to be preprocessed. 
Wiegelmann et al.l (120061 ) proposed a preprocessing routine to remove the net magnetic force 
and net magnetic torque, and to smoot h the magneti c field , while keeping the changes of 
the magnetic field within the noise level. iTadesse et al.l (|2009[ ) furt her developed t he pre pro- 
cessing routine into the spherical geometry. Following the idea of iTadesse et al.l (|2009[ ). we 
have written a version of the preprocessing code and applied it to the following analysis. 

We add Gaussian-distributed random noise to the bottom boundary of the Low and 
Lou solution to simulate realistic observations. For each component of the magnetic field 
{Br, Bq, or B^) at position i, a noisy component 6B^ = ■ {'-f = r,6 or 0) is added to the 
original magnetic field, where is the level of noise for that component, and rj is the random 
number in the range of —1 to 1 for position i. For all the grid points, r, has a Gaussian 
distribution. By choosing the noise level we have two different test cases here: 
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Noisy model 1: = 0.1 x A/niax(|i?,.|), Ig = ^ym.a,x{\Bg\), and = ^ymax{\B^\). 
Noisy model 2: Ir = 0.2 x A/max(|i?r|), = 2 x A/max(|i?e|), and 1^ = 2 x ^Jma^{\B^. 



The two noisy models are similar to noisy model 11 of IWiegelmann et al.l (120061 ). i.e., adding 
additional noises independent of the magnetic field strength. We only choose a different 
strategy to determine the noise level for each component of the magnetic field. The difference 
between noisy model 1 and noisy model 2 is that the latter has a noise level that is twice as 
large as that of the former. 

For each of the two noisy models, we apply the preprocessing procedure to remove 
the net magnetic force and net magnetic torque, and to smooth the magnetic field. The 
preprocessed bottom boundaries are submitted to the optimization procedure for NLFFF 
extrapolations. The other settings are the same as the test case 3 with an ideal bottom 
boundary, i.e., the initial condition and the boundary condition for the other five boundaries 
(other than the bottom one) are specified by the PFSS model. The extrapolation is also 
performed in the wedged shaped domain bounded by r G [1.0/?©, 2. Si?©], 6 G [9.0°, 171.0°], 
and (f) G [90.0°, 270.0°], which is resolved by 30 x 54 x 61 grid points. For comparison, 
the two noisy boundaries without preprocessing are also submitted for additional NLFFF 
extrapolations. The metrics for the extrapolation results are shown in the last four rows 
in Table [1] It shows that all the metrics for both noisy models 1 and 2 are improved after 
preprocessing of the bottom boundary. Since the normalized vector error metric and the 
mean vector error metric i?M are more sensitive to the difference of two sets of vector fields, 
we focus on the changes of these two metrics for the two noisy models before and after the 
preprocessing. The changes of E-^ and Em show that if the noises are larger in the bottom 
boundary, the preprocessing could improve the extrapolation results more significantly. 



3.4. Comparison between Cartesian and Spherical Extrapolations 

We cut a cubic box in the Low and Lou analytical magnetic field. The bottom surface 
of the box is tangent with the solar surface at [r, 6, 0] = [l.Oi?©, 78.0°, 180.0°]. The length of 
each side, /, is determined by I = 2i?Q tan(6'//2), where 9i is a free parameter indicating the 
angle of the box side viewed from the center of the Sun. For example. Figure [2] shows the case 
where 6i = 60°. We prepare three cases with 6i = 20°, 40°, and 60°. For each case, we do two 
NLFFF extrapolations in the spherical and Cartesian coordinates, respectively. The bottom 
boundary for the NLFFF extrapolation in the spherical geometry is the analytical magnetic 
field at IRq and bounded by the longitude and latitude lines at the four extremities of the 
bottom surface of the cubic box. The bottom boundary for the Cartesian extrapolation 
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is also the analytical magnetic field at l-R©, but the geometry of the magnetic fields are 
projected to the tangent bottom surface of the cubic box, while the vector components are 
kept in the heliographic coordinates. The other five boundaries for both cases are given by 
the potential field. Such a way to prescribe the boundary condition mimics real observations, 
since we only have real observations on the bottom boundary at present and people usually 
assume that the heliographic components are equivalent to the Cartesian components when 
the field of view is small. 

We adopt a buffer zone with 6 grid points for the weighting function. The preprocessing 
procedure is not applied to the boundaries. The NLFFF extrapolations are performed in 
the cubic box and the wedge shaped region for the Cartesian and spherical coordinate cases, 
respectively. The metrics for both the extrapolations in the Cartesian and spherical geome- 
tries are listed in Tabled The force-free measure, Lf, and the divergence-free measure, L^, 
are computed in the corresponding computation domain, i.e., the cubic box for the Cartesian 
case and the wedged shaped region for the spherical one. All the other metrics are computed 
in the inner region (excluding the six-grid buffer region) of the corresponding cubic box. 
Note that even for the spherical extrapolations, only the magnetic field in the inner region 
of the corresponding cubic box is extracted for computing the metrics. All the metrics in 
Table [2] show that the NLFFF code in the spherical geometry performs better than that in 
the Cartesian one when the field of view of the bottom boundary is larger than 20° x 20°. 

4. Applications of the Spherical NLFFF to Observations 

In order to apply the NLFFF model to observations, we need vector magnetic field 
measurements in a relatively large field of view and the initial condition that is often from 
the PFSS model, which is computed by the boundary condition from a global radial magnetic 
field. The S'DO/HMI provides us such magnetic fields with its full disk, high resolution, and 
high cadence vector magnetograms. 

4.1. Instrument and Data Analysis 

SDO /HMI observes the full Sun with a 4K x 4K CCD, whose spatial sampling is 0.5" per 
pixel. It obtains raw filtergrams at six different wavelengths and six polarization states (i.e., 
J ±5, where S* = Q, f/, and V) in the Fe I 6173 A spectral line. It takes 135 s to acquire each 
set of the filtergrams for deriving the vector magnetic field. All the filtergrams are averaged 
over 12 minutes to increase the signal to noise ratio and to remove the p-mode signals. The 
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four Stokes parameters, /, Q, U, and V, at the six wavelengths are computed from the raw 
data after necessary cahbrations, such as flat field correction, dark frame subtraction, and 
so on. 

The vector magnetic field, together with other thermodynamical parameters, are corn - 
puted by the code of Very Fast Inversion of the Stokes Vector (VFISV; iBorrero et al.ll201ll ). 
which fits the synthesized Stokes profiles to the observed ones by a least square fitting 
method. VFISV adopts a line synthesized model that is based on the Milne-Eddington at- 
mosphere, where all the physical parameters are constant along the line of sight, except for 
the source function, which depends linearly on the optical depth. VFISV has been optimized 
for the HMI data preprocessing; for instance, the damping constant has been set to 0.5 and 
the filling factor to 1. Therefore, the atmosphere is assumed to be uniformly filled with 
magnetic fields in each pixel. 

The transverse components of vector magnetic fields suffer from the so-called 180° ambi- 
guity. The 180° ambiguity for the HMI data in this study has been resolv e d by an improve d 



version of the m inimum ene r gy me thod (jMetcalllll994l : iMetcalf et al.ll2006l : iLeka et al. 



20091). 



As described in iLeka et al.l (120091 ). in weak-field areas, the minimization may not return a 
good solution due to large noise. Therefore, in order to get a spatially smooth solution in 
weak-field areas, we divide the magnetic field into two regions, i.e., strong- field region and 
weak-field region, which is defined to be where the field strength is below 200 G at the disk 
center, and 400 G on the limb. The values vary linearly with distance from the center to the 
limb. The ambiguity solution in the strong-field area is derived by the annealing in the min- 
imization method and released in the data series "hmi.B_720s_el5wl332_CEA" . Magnetic 
field azimuths in the weak-field finally determined by the potential-field acute-angle 

method after the annealing. This is different from the original minimum energy method, 
which uses a neighboring-pixel acute-angle algorithm to revisit the weak-field area. 

If a magnetic field vector is observed far from the solar disk center, the projection effect 
distorts the geometric shape of observed data, and the line-of-sight and the two transverse 
components {B^, Brj, and 5^) decomposed in the image plane coordinates deviate from the 
heliographic east-west, south-north, and vertical components {Bx,By, and B^). Therefore, 
correction for the projection effect must include two aspects. One corrects the geometric 



shape of the ob served data, which is often called remapping (ICalabretta fc GreisenI 12002 



Thompson! l2006l ) . The other transf orms the magnetic field components in the image plane 
to those in the heliographic plane (IGary fc HagyardI Il990l ). The heliographic components 
Bx,By, and B^ are equivalent to the spherical components B^, —Be, and Br. 

Figure [3] shows two sets of remapped and reprojected vector magnetic fields observed 
by HMI at 02:58 UT and 04:10 UT on 2012 January 23. The Joint Science Operations 
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Center (JSOC) series name for the data is "hmi.ME_720s_el5wl332" . More detail about 
the data are described in http://jsoc.stanford.edu/jsocwiki/VectorMagneticField An M8.7 
class flare peaked at 03:38 UT in active region NOAA 11402 on the same day. There are 
three additional active regions, NOAA 11401, 11405, and 11407 that interacted with NOAA 
11402 in this area. The field of view covering all the active regions is about 71.7° x 71.5°. 
Due to the large field of view, the NLFFF in the spherical geometry is necessary for studying 
the magnetic field configuration of the flaring active region and the magnetic connections 
between these active regions. 



4.2. Boundary and Initial Conditions 

Because vector-field measurements are not available for the side and top boundaries of 
a localized domain, one needs to make assumptions about these fields before performing a 
NLFFF extrapolation. Here, we assume the lateral and upper boundaries of the computa- 
tional domain are current-free, and make use of the well-known PFSS model to calculate the 
boundary conditions. In order to include as much information of magnetic field connections 
as possible and conform to the requirement of the "pfss" package, we need a global radial 
magnetic field as the bottom boundary for PFSS computations. But for realistic observations 
up to now, we could only get the full disk magnetogram, or the full disk vector magnetic 
field, in the best situation. It means that there are at least the following problems to be 
solved in order to construct a global radial magnetic field. 

First, we do not have magnetic field observations on the backside of the Sun. Therefore, 
the global magnetic field is usually constructed by the synoptic map, which is made by 
consecutive slices of magnetograms when each of them passes the central meridian. The 
slices of magnetograms are remapped on the heliographic plane, and are arranged side by side 
according to their Carrington coordinates. A synoptic map implicitly neglects the differential 
rotation of the solar surface, and only has a temporal resolution of the Carrington rotation 
period, i.e., approximately 27 days. In order to study the change of large-sca le magnetic 



structures on the time scale much less than the Carrington rotation period, IZhao et al. 



( 1l999l ) constructed a so-called "magnetic frame" consisting of two components, i.e., a classic 



synoptic map inserted with a remapped magnetogram observed at a short time range. 

Secondly, the magnetic field in the polar region of the Sun is not presently well observed. 
The poles are blocked periodically due to the tilt angle of about 7.25° between the Sun's 
equatorial plane and the ecliptic plane. The radial fields in the polar region are approxi- 
mately transverse to lines of sight as observed by instruments near the ecliptic plane, which 
introduces high noises. However, the large-scale polar magnetic fields are essential to build 
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the PFSS model; therefore, they are important for the NLFFF modehng. ISun et aLl (1201 if ) 
proposed a new two-dimensional spatial/temporal interpolation and flux transport model 
to interpolate the desired data. When we use the observed data (synoptic map or synoptic 
frame) as the boundary condition for the NLFFF extrapolation in the following sections, we 
adopt the new interpolation method to fill the magnetic field in the polar regions. 

Finally, a synoptic map or synoptic frame still lacks some basic physics of the magnetic 
field behavior on the solar surface, such as the differential rotation and the evolution of 



the rn agnetic field during a Carrington rotation. ISchrijverl (120011 ) and ISchrijver &: De Rosa 
( 120031 ) developed a flux- dispersal model and a data-assimilation procedure to address the 
aforementioned problems. The fiux-dispersal model advects the magnetic flux across the solar 
surface and incorporates all the ingredients that are necessary for the magnetic flux evolution, 
such as the differential rotation, meridional flow, convective dispersal, flux emergence, flux 
fragmentation, and flux cancellation. SDO /H.MI or the Michelson Doppler Imager (MDI) 
line-of-sight magnetograms are assimilated into the model when the data are available. 



We apply the PFSS model to two different kinds of photosp 



data as described above. One is the synoptic frame (IZhao et al. 



leric radial magnetic field 
I999I) a nd the other is an 



evolving flux dispersal model ( ISchrijverl I2OOII : ISchrijver fc De Rosa]l2003l ) that makes use of 
magnetograms assimilated into the model. These global magnetic maps are used to construct 
the PFSS models that are used for the initial flelds, and for the boundary conditions on the 
sides and top for both the selected observation times. The steps to construct the global 
radial magnetic fleld for each observation time are detailed as follows. 

First, the two vector magnetic fields as shown in Figure [3] are interpolated to 144 x 144 
grid points uniformly distributed in latitude and longitude with a grid spacing of about 0.5° 
per pixel. The vector magnetic fields are then preprocessed to remove the net magnetic force, 
torque, and noises using the preprocessing method in the spherical geometry as discussed in 
Section 13.31 

Secondly, the synoptic map for the Carrington Rotation 2119 (2012 January 9 to Febru- 
ary 6) is interpolated onto a uniform 720 x 360 pixel grid in latitude and longitude. The 
synoptic map is produced by line-of-sight magnetic field observations, which have been con- 
verted to the radial component by assuming that the magnetic field is radial. This synoptic 
map is used for both the observations at 02:58 UT and 04:10 UT, as we assume that the 
magnetic field outside the group of active regions did not change too much during this pe- 
riod. The synoptic map is aligned with the vector magnetic fields by a correlation method. 
Then, the region where the SDO /B.MI data are available is replaced by the preprocessed 
radial magnetic field. These new maps, or synoptic frames, are shown in Figures 111(a) and 
m^c) for the observations at 02:58 UT and 04:10 UT, respectively. 
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As an alternative, we prepare another version of the global radial magnetic field using an 
evolving flux dispersal model. For the observation at 02:58 UT, the surface radial magnetic 
field of the flux dispersal model at 00:04 UT is interpolated to 720 x 360 grid points, and 
aligned with the radial component of the preprocessed SDO /HMl vector magnetic field at 
02:58 UT. Additionally, we combine the aligned fiux dispersal field with the SDO /HMl 
radial field to include more accurate information where the observed data are available. 
Compared to the magnetic frame, the fiux dispersal model inserted with the observed data 
can be sampled every 6 hr and thus may provide a more consistent bottom boundary, since 
the fiux dispersal model has considered the differential rotation, meridional flow and other 
necessary factors. For the observation at 04:10 UT, we use the flux dispersal model at 
06:04 UT. Figures IH^b) andHt^d) show the flux dispersal model inserted with the SDO /HMl 
magnetogram for the observations at 02:58 UT and 04:10 UT, respectively. 

The PFSS models based on the two different boundary conditions at the two different 
times are computed with /max = 360. When the radius is large, a smaller /max is used. The 
correlation coefficients between the radial fields of the PFSS models (br) on the bottom 
boundary and the corresponding global radial fields (Br) for all the four cases are between 
0.93 to 0.96. The normalized errors of br and Br (defined as J2\^r — -Sr|/^|-Sr|) are 
between 0.27 to 0.34. Figure [5] shows a sample set of field lines traced through the PFSS 
models, which will later serve both as the side and top boundary conditions as well as the 
initial conditions for the NLFFF extrapolation in spherical geometry. Comparing the PFSS 
models at two observation times, we find that the PFSS models computed with the magnetic 
frame have evolved less than those computed with the fiux-dispersal model with the inserted 
magnetogram. The reason is that the synoptic map used for the two observation times is 
the same. Comparing the PFSS models with the two different boundaries at each time, we 
find that the polar fields computed with the magnetic frames are stronger than that with 
the other boundary, and this is why the "hairy Sun" images in Figure look different. The 
force-free, divergence-free measures, and the magnetic energy for the PFSS models are listed 
in Table [31 They have the same value for the two different boundaries for each observation 
time. 



4.3. Results 

The NLFFF models for both the observation times are derived via the optimization 
method with the boundary and initial conditions specified by the SDO /HMl vector magnetic 
fields and the PFSS models applied to both the synoptic frame and the fiux-dispersal model. 
Table [3] lists the force-free, divergence-free measures, and the magnetic energy contained in 
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these models. As discussed in Section \A.2\ using different boundaries for the PFSS model 
does not affect the force-free, divergence-free measures, and magnetic energy for the derived 
PFSS model. However, as we use these PFSS models for the NLFFF extrapolation, the 
flux- dispersal model gives better force-free, divergence-free measures, and larger magnetic 
energy than that computed from the synoptic frame for both the observation times. Since 
the potential field magnetic energy is the same, the free magnetic energy derived from the 
flux- dispersal model at 02:58 UT (2.0 x 10"^^ ergs) is thus larger than that from the synoptic 
frame (1.5 x 10'^^ ergs). Similarly, the magnetic energy at 04:10 UT from the fiux-dispersal 
model (1.2 x 10^^ ergs) is also larger than that from the synoptic frame (1.0 x 10^^ ergs). 
The difference in free energies between the two times, which we assume corresponds to the 
energy released by the M8.7 flare, is 0.8 x 10'^^ ergs or 0.5 x 10^^ ergs for the flux-dispersal 
model or the synoptic frame as the boundary conditions, respectively. 

The errors listed in Table |3] are estimated by a pseudo Monte Carlo experiment. For 
each of the four cases in Table [3l We add some Gaussian-distributed random noise to the 
global radial magnetic field. In order to focus on the effects of the two global magnetic field 
models, i.e., the synoptic map and the flux-dispersal model, the random errors are not added 
to the vector magnetic field. The error level ( the standard deviation of the random noises) 



is selected to be 10 G flHoeksema et al.ll2012l ). We create 10 noisy global radial magnetic 



fields for each case as listed in Table |3l combine them with the magnetogram observed by 
SDO /HMl, and do the PFSS and the NLFFF extrapolations. The metrics listed in Table [3] 
are the mean values derived from the extrapolation results computed from the 10 noisy and 
1 original boundaries. The errors are estimated as the standard deviations. 

We plot the PFSS and NLFFF models in Figure |6]in the wedge shaped field of view for 
the two observation times. Note that only the NLFFF extrapolations from the flux-dispersal 
model are plotted, since this model yields smaller force-free and divergence-free metrics. In 
these figures, the viewer is located above the center of the computational domain. For each 
observation time in Figure O the field lines of PFSS and NLFFF models are integrated from 
the same footpoints. We find that some field lines of the NLFFF at 02:58 UT before the M8.7 
flare expanded wider and higher than the PFSS model at the same time. A quantitative 
analysis is given below. Such a difference is not so obvious for the models at 04:10 UT. As 
for the magnetic connections between different active regions, we find that NOAA 11401 
plays a central role in this active region group. NOAA 11402 and NOAA 11401 form a long 
loop arcade. Some field lines emanating from NOAA 11401 fall into NOAA 11405. There 
are also some field lines emanating from the negative polarity of NOAA 11401 and 11407 
open to the interplanetary space. 

In order to evaluate whether the NLFFF models are consistent with the observed coronal 
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loops, and to gauge how much better the NLFFF models are when compared with the 
associated PFSS models, we plot the magnetic field lines over the SDO / AlA EUV images 
as shown in Figure [3 The center of the view point is located at (Lq, Bq), where Lq = 0.0° 
and Bq f» —5.3° are the longitude and latitude of the center of the solar disk, respectively. 
The composite EUV image from three AIA channels (211 A, 193 A, and 171 A, whose 
characteristic log(T) correspond to 6.3, 6.1, and 5.8, respectively) at 02:58 UT in Figure[71^a) 
shows that the magnetic loops in active region NOAA 11402 are highly nonpotential, since 
the distance between the footpoints are less than the height of the loops, and the core region 
is much brighter than the surroundings, which is an indication of large electric current. The 
PFSS field lines in Figure [7|(c) have an obvious deviation from the observed EUV loops, since 
the projection of the field lines leans towards the west, while the loops towards the east. The 
spatial correspondence between the overall shape of the NLFFF field lines and the EUV loops 
is much improved as shown in Figure [Tl^d). Therefore, the qualitative comparison between 
the model magnetic field lines and the observed EUV loops indicates that the NLFFF model 
provides a more consistent field for this active region. We recognize that the qualitative 
nature of this comparison is not ideal, however it remains the best option in the absence of 
a more reliable quantitative alternative. 

After the M8.7 fiare, the high EUV loops became invisible, while the post-fiare loop 
system was more conspicuous as shown in Figure [Tl^b). The difference between the magnetic 
field lines of the PFSS and NLFFF models at 04:10 UT is not so obvious as shown in 
Figures Wie) and (f), which is consistent with the fact that the electric current has been 
ejected along with the fiare and filament eruption. The NLFFF at 04:10 UT contains less 
free magnetic energy than that at 02:58 UT as listed in Tabled which provides additional 
evidence that the overall magnetic field configuration relaxed to a more potential state after 
the fiare. 

In order to give a quantitative comparison between the PFSS and NLFFF field lines, we 
measure the distance of field line footpoints and the maximum height of a field line. The ratio 
between them tells roughly the shape of a field line. The sample field lines are the same as 
those in Figures [6] and [71 but emanating from the active region NOAA 11402. The distance 
is measured along the great circle on the solar surface. Figure [H] shows the footpoint distance 
versus the maximum height of the field line for both the PFSS and NLFFF models before 
and after the M8.7 fiare. For the field lines higher than 50 Mm above the solar surface at 
02:58 UT, there are quite some NLFFF field lines whose maximum heights are greater than 
their footpoint distances, but there are very few PFSS field lines having this property. Such 
a behavior becomes weaker for the field lines at 04:10 UT. Therefore, we have a quantitative 
measurement that some NLFFF field lines in the fiare productive region expanded higher 
than the PFSS field lines before the fiare. 
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5. Summary and Discussions 

We test a NLFFF optimization code in the splierical geometry witli tlie Low and Lou 
analytical solution. It is found that the analytical solution can be recovered to a high level 
of accuracy if vector data on all boundaries are provided. For a more realistic simulation, 
we only provide the vector magnetic field on the bottom. The results show that a weight- 
ing function for the lateral and top boundaries and a pre-processing for the noisy bottom 
boundary are necessary for a reliable NLFFF extrapolation. Analytical tests also show that 
the NLFFF code in the spherical geometry performs better than that in the Cartesian one 
when the field of view of the bottom boundary is large, say, 20° x 20°. 

We apply the NLFFF model to an active region observed by SDO /HMl both before and 
after an M8.7 flare. The initial condition for the NLFFF optimization method in the spherical 
geometry is computed by the PFSS model. PFSS uses a global radial magnetic field as its 
boundary condition, which can be derived by different methods, such as the synoptic frame 
or the flux- dispersal model. For each observation time, we compare the two results computed 
with different initial conditions derived by the PFSS model using the synoptic frame and 
the fiux-dispersal model, respectively. The results show that NLFFF extrapolations using 
the fiux-dispersal model as the boundary condition (to compute their initial conditions) have 
slightly lower, therefore better, force-free and divergence-free metrics, and contain larger free 
magnetic energy. However, if we consider the errors as listed in Table E] and discussed in 
Section 14. 3[ the metrics are marginally comparable. 

We compare the extrapolated magnetic field lines both from the NLFFF and PFSS 
models with the EUV observations by SDO / AlA. The comparison shows that NLFFF per- 
forms better than PFSS not only for the core field of the flare productive region but also for 
the large EUV loops higher than 50 Mm. On the one hand, NLFFF bears larger magnetic 
energy than PFSS. On the other hand, the EUV loops in active region NOAA 11042, where 
the M8.7 flare occurred, have non-potential shapes. NLFFF field lines reproduce these EUV 
loops better than PFSS before the flare. In this active region group, we have both closed 
and open field lines, both bipolar and quadrupolar magnetic configurations. Thus, it is a 
good example for us to further study the flare itself and its associated filament eruption and 
CME. After the flare, the difference between NLFFF and PFSS becomes smaller than that 
before the flare. 

Compared to previous studies on NLF FF extrapolations in spherical geometry (e.g.. 



Wiegelmannl 120071 : iTadesse et al.ll2009l . 120111 ). the following additional experiments are con- 
ducted in this study. First, we quantitatively compare the extrapolation results derived by 
the NLFFF codes both in the spherical and Cartesian coordinates using the Low and Lou 
analytical solution. These comparisons show that geometrical factors need to be taken into 
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account when the field of view of the bottom boundary is large. Next, we have initialized 
the NLFFF calculations with global PFSS extrapolations, in light of the fact that the active 
regions are rarely isolated and the global PFSS extrapolations provide some information 
about the connectivity of field lines between the active region of interest and points outside 
of the computational domains considered here. Additionally, we test the sensitivity of the 
resulting NLFFF models on the surface fields, and perform comparisons between NLFFF 
extrapolations based on different surface fields (the synoptic frame and fiux-dispersal model). 
Finally, we apply the NLFFF extrapolation in spherical geometry to SDO /U.M1 observations 
and compare the magnetic field lines with EUV coronal loops observed by SDO / AlA. We 
conclude that the NLFFF model performs better than the PFSS model, especially for highly 
non-potential active regions before eruptive activities. 

Due to the limitation of computation resources, we adopt a spatial resolution of 0.5° 
(note that 0.5° is the longitude or latitude difference) per grid point for this study. Since 
we focus on the overall property of the core field region, such as the magnetic energy and 
the large scale magnetic loops, such a resolution is adequate. However, the SDO/BMl has 
a much higher resolution, which is approximately equivalent to 0.03° per pixel. In future 
research, if one needs to explore the small-scale structure, it would be better to interpolate 
the initial condition to a higher resolution, combine it with the SDO /HMl vector magnetic 
fields, and then perform the NLFFF extrapolation to increase the spatial resolution. 



Recently, IWiegelmann fc Inhesterl (120101 ) added a new term in the opt imization method 



to dea l with measurement errors and missing data on the bottom boundary. IWiegelmann et al 



(120121 ) ha s tested this op t imiza tion method in the Cartesian geometry with the SDO observa- 



tions, and lTadesse et al.l ( 12012! ) applied the spherical version to observations from the Synop- 
tic Optical Long-term Investigations of the Sun (SOLIS). The new optimization method has 
the advantage to deal with measurement errors and minimize the force-free and divergence- 
free terms, while it takes much longer time for computation than the traditional optimization 
method. Although the NLFFF extrapolation techniques have been improved in various ways, 
it is still a challenge to validate all these magnetic field models with coronal observations. 
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Table 1: Force-free, divergence- free, vector correlation, Cauchy-Schwartz, normalized vector 
error, and mean vector error metrics and normalized magnetic energy for the Low and Lou 
solution, potential field source surface model, and nonlinear force-free field models. 







T 2 


'-'vcc 


2 


1 




1 — Eyi^ 


e2 


Model 1 


(10-2 


G2 Mm-2) 














Low & Lou (1990) 


0.03 


0.01 


1 


1 


1 




1 


1 


PFSS (/max = 7) 


0.01 


0.01 


0.85 


0.83 


0. 


51 


0.45 


0.82 


Case 1 


0.00 


0.00 


1.00 


1.00 


0. 


99 


0.99 


1.01 


Case 2 


0.34 


0.11 


0.98 


0.90 


0. 


78 


0.61 


0.92 


Case 3 


0.01 


0.00 


0.99 


0.95 


0. 


86 


0.72 


1.01 


Noisy model 1 


0.24 


0.16 


0.99 


0.94 


0. 


81 


0.68 


0.97 


Preprocessed 1 


0.04 


0.02 


0.99 


0.95 


0. 


85 


0.72 


1.00 


Noisy model 2 


0.91 


0.61 


0.97 


0.90 


0. 


72 


0.59 


0.99 


Preprocessed 2 


0.11 


0.06 


0.99 


0.94 


0. 


84 


0.71 


1.00 


Notes. 



^ See Section 13.31 for the definitions of the models. 

2 Refer to equations (l6|)- (fT2l) in Section 13.31 for the definitions of the metrics. All 
the metrics are computed in the computation domain of r G [l.Oi?©, 2.5i?0], 9 G 
[9.0°, 171.0°], and G [90.0°, 270.0°]. 
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Table 2: Force-free, divergence- free, vector correlation, Cauchy-Schwartz, normalized vector 
error, and mean vector error metrics and normalized magnetic energy for the nonlinear 
force-free field models in Cartesian and Spherical Geometries. 









La' 


r 3 

'-'vcc 




1 — 


1 — Eu^ 




Model 1 




(10^2 


Mm-2) 












Cartesian ( 


[01 = 20°) 


96.92 


57.91 


0.80 


0.65 


0.42 


0.33 


0.68 


Spherical ( 


ei = 20°) 


0.36 


0.21 


0.95 


0.89 


0.71 


0.63 


1.06 


Cartesian ( 


% = 40°) 


21.98 


18.02 


0.89 


0.80 


0.50 


0.40 


0.98 


Spherical ( 


9i = 40°) 


0.16 


0.12 


0.98 


0.94 


0.82 


0.72 


0.99 


Cartesian I 


% = 60°) 


6.62 


5.80 


0.90 


0.80 


0.42 


0.37 


1.64 


Spherical ( 


01 = 60°) 


0.09 


0.07 


0.98 


0.92 


0.79 


0.65 


0.93 



Notes. 

^ See Section 13.41 for the definitions of the models. 



^ Refer to equations ([6]) and ([7]) in Section 13.31 for the definitions of the metrics. Lf 
and Ld are computed in the corresponding computation domain. 
3 Refer to equations fl8|)-f lT2|) in Section 13.31 for the definitions of the metrics. These 
metrics are computed in the inner region (excluding a six-grid buffer region) of the 
corresponding cubic box extracted from their computation domain. 



Table 3: Force-free and divergence-free metrics and magnetic energy for the nonlinear force-free field models and 
potential field with observed data as the boundary condition. 





Time 




L^ 




2^poti 


rpoti 


pnlfff2 


ppot2 
m 


77'free2 


Date 


(UT) 


Boundary 


(10-2 G2 Mm-2) 






(10^3 ergs) 






2012 Jan 23 


02:58 


synoptic frame 


3.53 ±0.06 


1.78 ±0.01 


0.25 ±0.00 


0.26 ±0.00 


1.58 ±0.02 


1.43 ±0.00 


0.15 ±0.02 


2012 Jan 23 


02:58 


flux-dispersaP 


3.41 ±0.08 


1.76 ±0.02 


0.25 ±0.00 


0.26 ±0.00 


1.63 ±0.03 


1.43 ±0.00 


0.20 ±0.03 


2012 Jan 23 


04:10 


synoptic frame 


4.76 ±0.10 


2.20 ±0.02 


0.32 ±0.00 


0.33 ±0.00 


1.67 ±0.01 


1.57 ±0.00 


0.10 ±0.01 


2012 Jan 23 


04:10 


fiux-dispersaP 


4.65 ±0.13 


2.18 ±0.03 


0.32 ±0.00 


0.33 ±0.00 


1.69 ±0.01 


1.57 ±0.00 


q1^2±0.01 

— CT3 



Notes. 



Lf, Ld, L\° , and are the force-free, divergence-free measures for the nonlinear 
force-free field and the potential field, respectively. These metrics are computed in 
the computation domain including the buffer region. 

2 -^m°*5 E^^'^ are the nonlinear force-free field, potential field, and the 
free magnetic energy, respectively, where E^^^ = E^^^ — E^^. These metrics are 
computed in the inner region (excluding a six-grid buffer region) of their computation 
domains. 

3 The observed magnetic field has been inserted to the fiux-dispersal model. 
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Fig. 1. — (a) Low and Lou nonlinear force- free field solution. Red, green, and black lines 
represent open field lines rooted on negative polarities, open field lines rooted on positive 
polarities, and closed field lines, respectively. The blue line indicates a wedge shaped domain 
surrounded by r e [1.0Rq,2.5Rq], 9 e [9°, 171°], and G [90°, 270°]. All the metrics in 
Table [1] are computed in this domain, (b) Potential field source surface model computed 
with the boundary of the Low and Lou solution, (c) Nonlinear force-free field computed 
with the noisy boundary 1 (see texts for details), (d) Nonlinear force- free field computed 
with the preprocessed noisy boundary 2. 
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Fig. 2. — A wedge shaped region cut from the Low and Lou solution. The center of the 
bottom surface is located at [9, cf)] = [78°, 180°]. Blue lines mark the boundary of the wadge 
shaped region. Cyan lines represent a cubic box that is tangent to the solar surface at 
[r,9,(j)] = [li?0, 78°, 180°]. Red, green, and black lines represent open field lines rooted 
on negative polarities, open field lines rooted on positive polarities, and closed field lines, 
respectively, (a) The view angle is [^,0] = [100°, 170°]. (b) The view angle is [6*, 0] = 
[70°, 90°]. 
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SDO/HMI 6173 23-Jan-2012 02:58:10.700 UT SDO/HMI 6173 23-Jan-2012 02:58:10.700 UT 
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Longitude (degree) Longitude (degree) 



Fig. 3. — Vector magnetic fields observed by SDO /B.MI. (a) The vertical component of 
the vector magnetic field observed at 02:58 UT on 2012 January 23. The black rectangle 
represents the field of view for panel (b). Written in the rectangle are the NOAA numbers 
for the active regions, (b) The vector magnetic field observed at 02:58 UT in the field of view 
as shown in the rectangle of panel (a). The arrows stand for the horizontal components, (c) 
The vertical component of the vector magnetic field observed at 04:10 UT. (d) The vector 
magnetic field observed at 04:10 UT. 
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Synoptic Frame 23-Jan-2012 02:58:10.700 UT Flux Dispersal + Magnetogram 23-Jan-2012 02:58:10.700 UT 



I 



(b) 



1 



-100 100 -100 100 

Longitude (degree) Longitude (degree) 

Synoptic Frame 23-Jan-2012 04:10:10.700 UT Flux Dispersal + Magnetoqram 23-Jan-2012 04:10:10.700 UT 





Longitude (degree) 



-100 
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Fig. 4. — (a) The synoptic frame that is a combination of the synoptic map for the Car- 
rington Rotation 2119 (2012 January 9 to February 6) and the magnetogram observed by 
SDO/BMl at 02:58 UT on 2012 January 23. The black rectangle indicates the field of view 
of the magnetogram. (b) The global vertical magnetic field that is a combination of the 
flux dispersal model at 00:04 UT on 2012 January 23 and the magnetogram observed by 
SDO /HMl at 02:58 UT on 2012 January 23. (c) Similar to panel (a) but using the magne- 
togram at 04:10 UT. (d) Similar to panel (b) but using the flux dispersal model at 06:04 UT 
and the magnetogram at 04:10 UT. 



- 31 - 




Fig. 5. — The PFSS models computed with the four different boundaries as shown in FigureHl 
Panels (a)-(d) and those in Figure H] have a one-to-one correspondence. Red, green, and 
black hues represent open field lines rooted on negative polarities, open field lines rooted on 
positive polarities, and closed field hues, respectively. 
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Fig. 6. — (a) The PFSS model extracted from panel (b) of Figure [5] in the domain of 
r e [1.0Rq,2.5Rq\, 9 E [29.1°, 100.8°], and E [-5.0°, 66.5°]. The yellow lines indicate 
the magnetic field lines, whose maximum heights are less than l.O6i?0. The cyan lines are 
for other closed magnetic field lines. The red lines represent open magnetic field lines, (b) 
The NLFFF model computed by the vector magnetic field observed at 02:58 UT and the 
PFSS model in panel (a), (c) The PFSS model extracted from panel (d) of Figure O in the 
domain of r G [l.Oi?©, 2.5/?©], 9 E [29.1°, 100.8°], and E [-4.5°, 67.0°]. (d) The NLFFF 
model computed by the vector magnetic field observed at 04:10 UT and the PFSS model in 
panel (c). 
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Fig. 7. — (a) The EUV image at 211 A , 193 A , and 171 A bands (corresponding to red, 
green, and blue channels, respectively) observed by SDO/AIA at 02:58 UT on 2012 January 
23. (b) The EUV image observed at 04:10 UT. (c) Field lines of the PFSS model at 02:58 
UT overlaid on the EUV image. The color code is the same as that in Figure El (d) Field 
lines of the NLFFF model at 02:58 UT overlaid on the EUV image, (e) Field lines of the 
PFSS model at 04:10 UT overlaid on the EUV image, (d) Field lines of the NLFFF model 
at 04:10 UT overlaid on the EUV image. 
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Distance (Mm) Distance (Mnn) 



Fig. 8. — The distance of the two footpoints versus the maximum height of a magnetic field 
hne. The asterisks are for the case of the PFSS field lines, and the squares for the NLFFF field 
lines. The dash-dotted line marks the locations where distances equal to heights, (a) The 
magnetic field lines are computed with the data at 02:58 UT, and integrated from the same 
footpoints as shown in Figures Et^a), Mh), C^c) and [7(d), but in the region 9 G [50.0°, 70.0°] 
and (p G [10.0°, 30.0°]. (b) The magnetic field lines are computed with the data at 04:10 UT, 
and integrated from the same footpoints as shown in Figures [6t^c), [6](d), [7](e) and[7](f), but 
in the region 9 e [50.0°, 70.0°] and (p E [10.5°, 30.5°]. 



